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EVOLUTION OF MAGNETIC FIELDS IN SUPERNOVA REMNANTS 

K.M. Schure 1 , J. Vink \ A. Achterberg 1 and R. Keppens x ' 2 ' 3 

RESUMEN 

En este momento es generalmente aceptado que los remanentes de supernova (SNR) son la fuente de los 
rayos cosmicos (CRs) hasta una energia de 10 15 eV. Los campos magneticos requeridos para acelerar los CRs 
a energias suficientemente altas necesitan ser mucho mayores de los que pueden resultar de la compresion del 
medio circunestelar (CSM) por un factor 4, como ocurre choques fuertes. Mapas de emision sincrotron de estas 
regiones indican que de hecho el campo magnet ico es mucho mayor, y que tiene una configuracion radial para 
remanentes jovenes y toroidal para remanentes viejos. Tanto el mecanismo de fortalecimiento del campo como 
la razon de su configuracion radial (para remanentes jovenes) no son entendidos en forma satisfactoria. Usamos 
el codigo MHD de malla adaptativa AMRVAC para simular la evolucion de remanentes de supernova, y ver 
si podemos reproducir la configuracion de campo magnetico radial en los estados evolucionarios tempranos. 
Seguimos la evolucion del SNR con tres configuraciones iniciales del campo magnetico en el CSM: un campo 
inicialmente toroidal, un campo turbulento, y un campo paralelo al eje de simetria. A pesar de que en los 
ultimos dos casos una componente radial apreciable se genera en la discontinuidad de contacto como resultado 
de la inestabilidad de Ray leigh- Taylor, no se genera una componente de campo radial en la zona inmediatamente 
interior al choque externo. Aparentemente, la MHD ideal no parece producir el efecto observado. Posiblemente, 
una mayor compresion y turbulencia adicional debida a la presencia dominante de CRs nos podrian ayudar a 
reproducir mejor las observaciones. 

ABSTRACT 

Supernova remnants (SNR) are now widely believed to be a source of cosmic rays (CRs) up to an energy of 
10 15 eV. The magnetic fields required to accelerate CRs to sufficiently high energies need to be much higher than 
can result from compression of the circumstellar medium (CSM) by a factor 4, as is the case in strong shocks. 
Non-thermal synchrotron maps of these regions indicate that indeed the magnetic field is much stronger, and 
for young SNRs has a dominant radial component while for old SNRs it is mainly toroidal. How these magnetic 
fields get enhanced, or why the field orientation is mainly radial for young remnants, is not yet fully understood. 
We use an adaptive mesh refinement MHD code, AMRVAC, to simulate the evolution of supernova remnants and 
to see if we can reproduce a mainly radial magnetic field in early stages of evolution. We follow the evolution 
of the SNR with three different configurations of the initial magnetic field in the CSM: an initially mainly 
toroidal field, a turbulent magnetic field, and a field parallel to the symmetry axis. Although for the latter two 
topologies a significant radial field component arises at the contact discontinuity due to the Ray leigh- Taylor 
instability, no radial component can be seen out to the forward shock. Ideal MHD appears not sufficient to 
explain observations. Possibly a higher compression ratio and additional turbulence due to dominant presence 
of CRs can help us to better reproduce the observations in future studies. 
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1. INTRODUCTION 

Supernova remnants (SNRs) are interesting ob- 
jects in which particles are accelerated and vari- 
ous magnetic-field morphologies have been detected. 
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The widths of the observed X-ray synchrotron fil- 
aments indicate that the magnetic field is much 
stronger than can be expected from compression of 
the interst ellar field by the supernova blast wave 
alone (e.g. IVolk et alJbood IVinkl l2QQ5h . and the 



spectral properties indicate that the magnetic field 
turbulence is high enough that cosmic ray diffu- 
sion must be near optimal for cosmic ray a c celera - 
tion (the so - called Bohm-diffusion; e.g. lVinkl ( 2005 ): 



Stage et al.l (|2006f )). Polarization measurements in 



radio show that young supernova remnants have a 
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predominantly radial magnetic field, whereas the 
dominant orientation of the f ield in old er remnants 
is parallel to the shock front (|Gulllll973h . 

The radial field in young SNRs is thought to arise 
because of the Rayleigh-Taylor (R-T) instability that 
occurs at the contact discontinuity, which separates 
the shocked ejecta from the shocked, more tenu- 
ous circumstellar medium (CSM). The R-T insta- 
bility creates fingers into the more tenuo us medium , 
which stretch the magnetic field radially (| Gullll 1 9 73h . 
Magnet oh y dro dy namic (MHD) simulations by e.g. 
Jun et all ( 19951 ) have indeed shown that the R-T 
fingers cause a preferred polarization of the field 
along the fingers. While this causes the dominance of 
the radial field signature in part of the remnant, for 
a typical shock compression ratio of 4, it does not ex- 
plain the radial field polarization observed near the 
forward shock. As shown i n hydr odynamic simula- 
tions by iBlondin fc Ellisonl ( 200lf h when the com- 
pression ratio is higher, the R-T fingers can almost 
reach, and perturb the forward shock. The presence 
of cosmic rays could justify an adiabatic index lower 
than that for a monatomic gas, where 7 = 5/3, and 
thus a higher compression ratio (s = (7+ 1)). 
For a gas that is dominated by relativist ic particles, 
such as cosmic rays, 7 = 4/3. Additionally, the es- 
cape of high energy cosmic rays (CRs) drains energy 
from the shock region, thus further lowering the adi- 
abatic index. The signature toroidal magnetic field 
of older SNRs is thought to arise by compression of 
the circumstellar magnetic field. 

We explore the development of Rayleigh-Taylor 
fingers and magnetic field variations of a SNR inside 
a stellar wind medium (p oc r -2 ) for different initial 
magnetic topologies, and compare our results with 
those from various studies that have been performed 
on this subject in a homogeneous and/or unmag- 
netized medium. One of the magnetic field topolo- 
gies we will consider is that of a magnetized stellar 
wind. Even for a slowly rotating star, the field will 
be mainly toroidal at a distance much larger than the 
stellar radius. Another topology we consider is one 
of magnetic turbulence. The interstellar magnetic 
field has a large turbulent component, so this is a 
very relevant scenario for a SNR in a homogeneous 
medium. In a CSM created by a stellar wind how- 
ever, a large turbulent component could be driven by 
cosmic ray streaming ahead of the shock front and 
small perturbations of the blast wave. Alternatively, 
an ordered magnetic field, as observed on galactic 
scales, can be used to initialize the magnetic topol- 
ogy. This will be the third case we consider. For 
now, we focus on the early stages of evolution and 



try to see whether or not a radial field component is 
dominant in regions of a young SNR. 

2. METHODS 

We model the evolution of a SNR in a CSM that 
is shaped by a pre-supernova wind, in this case a 
slow, cool wind, such as that originating from a Red 
Supergiant (RSG), before explosion. The adopted 
mass loss rate is M = 1.54 x 10 -5 M yr _1 , the 
wind velocity v = 4.7 km s _1 and the temperature 
T = 1000 K. We use the Adaptive Mesh Refinement 
version of the Versatile Advect ion Code: AMRVAC 
van der Hoist fc Keppensll2007h to solve the MHD 



equations in the r — plane of a spherical grid, with 
symmetry around the polar axis. The grid spans 
1.8 x 10 19 cm radially, and the full 7r rad angle. Since 
the supernova remnant does not reach beyond the ra- 
dius where the red supergiant (RSG) wind meets the 
main sequence bubble in the timescales we consider, 
the initial grid is filled with a RSG wind. This is 
done by inducing the stellar wind at the inner radial 
boundary, and let this evolve until it fills the grid 
and a stationary situation has been established. We 
then artificially introduce the magnetic field in the 
entire grid, and supernova ejecta into the inner 0.3 
pc. 

The initial ejecta density profile consists of a 
constant-density core, with an envelope for which 
the density decreases as p oc r -9 , which is the 
typical density profile of the ejecta after propa- 
gation through the star in explosion models (c.f. 
iTruelove In order to match the 

observationally determined ejecta mass and energy, 
in this case chosen to be M e j = 2.5M and E e - } = 
2 x 10 51 erg, we iteratively determine the value for 
the density in the ejecta core and the velocity at 
which the core ends and the powerlaw envelope be- 
gins. For the parameters adopted in our simulations, 
this happens for a central density of the ejecta core of 
5.1 x 10 _20 g cm -3 , and the powerlaw envelope begins 
where the velocity of the ejecta is 8.7 x 10 8 cm s _1 . 
The velocity linearly increases from zero at the core, 
to 15,000 km s _1 at the outer part of the envelope. 

The number of gridzones in the simulations is 
300x90 on the first level, with 3 refinement lev- 
els (with ratio 2) for the wind-stage, and 5 refine- 
ment levels (refinement ratio 2) for the SNR stage. 
The refinement is based on the density and velocity 
of the fluid. We thus reach an effective resolution 
of 4800x1440, corresponding to 2.75 x 10 15 cm by 
0.125°, in the regions where strong density and ve- 
locity gradients are present. We implement three 
toy models with different magnetic field topologies 
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in the wind and track what happens when the SNR 
evolves. In all cases the ratio of magnetic energy 
density to kinetic energy density is about 0.01. This 
is low enough that the equilibrium wind solution is 
not significantly perturbed on time-scales relevant 
for our supernova remnant evolution. 

In Model A, we set up the field as if the RSG 
wind were magnetized and the star slowly rotating. 
The field at r ^> R s tar is predomina ntly toroidal, but 
has a very small radial componen t ( Chevalier fc Luol 



1994; iGarcfa-Segura et al 



ponent ([< 
0993). i 



-l 



e. B<f) oc r 
the toroidal field 



B r oc r . Far from the star 
dominates. Since the toroidal field has to vanish at 
the poles, w e induce a sinfl dependenc e, similar to 
that used bvlGarcfa-Segura et al. (119991) . It was sug- 
gested by iBiermann fc Cassinellil (|l993h that such a 
progenitor magnetic field may be responsible for high 
magnetic fields in SNRs, leading to efficient cosmic 
ray acceleration. 

In Model B, we set up a random 2D magnetic 
field in the r and directions at the time when we 
introduce the eject a. The component is kept 
zero. At all times, we implement an initial analyt- 
ically divergence- free magnetic field. The random 
fiel d is calculated on a 2D ca rtesian grid, follow- 
ing [Giacalone^^Joki^ (|l999l ). and transcribed to 
spherical coordinates in the r, plane. In a centered 
difference evaluation, we find that the initial condi- 
tion of the magnetic field is close to divergence-free 
(|(V-£)|/|B| <C l/\v£d\). The field is set up accord- 
ing to: 



5B(x,y) 



0.5a 



i(cos (j) n y - sin </> n x) 

ik n (x cos <\>n J rV si* 1 0n)+*/^n 



(1) 



with 256 different wavenumbers /c, and the phase 
(3 and polarization <p are randomly chosen between 
and 2tt for each wavenumber. The smallest 
wavenumber spans the entire grid. The wavenum- 
bers are logarithmically spaced, and the largest 
wavenumber represents a wavelength that covers two 
gridcells at the lowest resolution. For an approxima- 
tion of a Kolmogorov spectrum, we set a = 5/3. The 
above equation gives us a complex value for the mag- 
netic field. To get a real magnetic field, we take the 
real part of the above equation, by tacitely assum- 
ing that the complex part of Ao has been combined 
with the phase in (3. A turbulent magnetic field may 
be a natural outcome of cosmic ray induced mag- 
netic field amplification. Since the amplification is 
stronger around the fast shocks of young SNRs, this 
may help explaining why only young SNRs have ra- 
dial magnetic fields, if turbulent magnetic fields in 



front of the shock lead to radial magnetic fields in- 
side the SNR. 

In Model C, the magnetic field is set up paral- 
lel to the rotation axis, representing a situation like 
the dominant ordered magnetic field observed, e. g., 
in spiral galaxies. In spherical coordinates, this is 
represented by B r = BcosQ, Bq = —BsinO. 

In all simulations, the divergence of the mag- 
netic field is controlled by adding a source term pro- 
portional to V • B in the induction equation, while 
maintaining conservation of momentu m and energy 
([Keppens et al.ll2003l : IJanhunenllioOoL 



3. RESULTS 

The results from simulations of Models A, B and 
C, are plotted in Figures [T] to 03 The density, radial 
and toroidal magnetic field, and additionally the ra- 
tio of the radial to the total field are plotted for the 
three magnetic field topologies. The results shown 
here are for remnants that have evolved for a period 
of 634 year. The radial scale is 10 18 cm. 
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Fig. 1. Simulation of SNR into CSM with a mostly 
toroidal initial field (Model A). The upper left quadrant 
shows the logarithm of the density, which shows clearly 
the occurence of the Ray leigh- Taylor instability at the 
contact discontinuity. The lower left panel shows the 
logarithm of the absolute value of the radial magnetic 
field. The lower right part shows the logarithm of |i^|, 
and the upper right quadrant shows the logarithm of the 
ratio of the radial field component relative to the total 
Beld (|B r |/|B|). 
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Fig. 2. Simulation of the SNR into CSM with a turbulent 
initial field (Model B). Similar to Figure [1] the density 
and magnetic field is plotted. Note that in this case the 
lower right part shows the logarithm of \Be\, and the 
upper right quadrant shows the ratio of the radial field 
component relative to the total field (|5 r |/|B|). 

The supernova eject a sweeps up the circumstellar 
matter into a dense shell. Four distinct regions can 
be identified: the unshocked circumstellar medium 
(CSM) ahead of the blast wave, the shocked CSM, 
the shocked eject a, and the freely expanding eject a. 
The shocked CSM is separated from the unshocked 
CSM by a strong shock (in Figures [T]|3] at a radius 
of Rf « 1.3 x 10 19 cm), characterized by a pres- 
sure jump, which increases the density by a factor 
(7+I)/ (7 — 1) and decreases the velocity in the shock 
frame by the same amount. The shocked CSM is sep- 
arated from the third region, the shocked ejecta, by 
the contact discontinuity, characterized by a jump 
in the density but constant pressure. This is also 
where the R-T instabilities arise. The reverse shock, 
at a radius of R r w 6.5 x 10 18 cm, marks the bound- 
ary with the unshocked ejecta. The deceleration of 
the shocked ejecta by the less dense shocked CSM is 
Ray leigh- Taylor unstable. Since the magnetic field 
is very weak, it does not influence the dynamics of 
the blast wave, and we do not see a difference in the 
development of the R-T instabilities and the propa- 
gation of the blast wave between the three models. 

The magnetic field is carried along by the plasma. 
Since the velocity field initially only has a radial 
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Fig. 3. Simulation of SNR into CSM with an initial mag- 
netic field parallel to the symmetry axis, i.e. vertically 
aligned (Model C). The upper quadrants show again the 
logarithm of the density and the ratio of radial field over 
the total field. The lower quadrants show the absolute 
values of B r and Bq. 

component, the induction equation gives us that the 
toroidal field is swept up, while the radial field re- 
mains: <9tB = V x (v x B). Not surprisingly there- 
fore, in the ejecta, the radial field component is the 
dominant one. Starting at the R-T region at the 
contact discontinuity, the toroidal field starts to be 
present again. The R-T instability induces a vq com- 
ponent, which subsequently causes changes also in 
the radial component of the magnetic field. Along 
the R-T fingers, for Models B and C, the field is 
mostly radial with toroidal components at the tips 
and base s of the fingers, in agreement with results 
from e.g. Uun et al.l ([1995). In Model A however, we 
do not see a significant radial component in the R-T 
unstable region, which may have to do with the small 
scale of the initial radial field, something requiring 
further investigation. 

The radial field along the R-T fingers in Models B 
and C is not sufficient to explain the observations of 
a dominant radial polarization of the field in young 
SNRs up to the shock front. In Model B, a few fea- 
tures in which the field is predominantly radial can 
be seen out to the forward shock, but not enough to 
explain observations. As the R-T is not yet fully sat- 
urated at this time in the evolution, the radial field 
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may stretch out closer to the forward shock. How- 
ever, since in very young remnants the field is al- 
ready predominantly radial, this alone is insufficient 
to explain observations. 

4. CONCLUSIONS AND DISCUSSION 

Unlike what is observed, the radial field does not 
dominate in the outer region of the remnant in our 
models. Although in model B the field fluctuates 
more due to the random initial condition, there is 
no distinct radial field in the remnant. Specifically 
with an initial toroidal field, it is difficult to create 
a field that is mainly radial. Turbulence could be an 
essential component in creating a radial field, but as 
shown in the result of model B, this is by itself not 
sufficient to reproduce the observed geometry of the 
magnetic field of young remnants. 

It appears that additional physics is needed to 
explain the observations of a high radial field in 
young SNRs, such as a softened equati on of state as 
suggested by iBlondin fc Ellison! (|200ll ) : a situation 
that may occur when cosmic rays become dynami- 
cally important at the shock. A turbulent magnetic 
field may arise and be amplified by CR streaming 
(|Bell fc Lucekl lioQlh . The turbulence in turn, can 
increase the maximum energy that CRs can attain. 
In future studies we plan to explore the interaction 
between cosmic rays and magnetic fields in super- 
nova remnants. 
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